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Greatly Enhanced Eccentricity Oscillations in Quadruple Systems 
Composed of Two Binaries: Implications for Stars, Planets, and 
Transients 
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ABSTRACT 

We study the orbital evolution of hierarchical quadruple systems composed of two binaries on 
a long mutual orbit, where each binary acts as a Kozai-Lidov (KL) perturber on the other We 
find that the coupling between the two binaries qualitatively changes the behavior of their KL 
cycles. The binaries can experience coherent eccentricity oscillations as well as excursions to 
very high eccentricity that occur over a much larger fraction of the parameter space than in 
triple systems. For a ratio of outer to inner semi-major axes of 10 to 20, about 30 to 50% of 
equal-mass quadruples reach eccentricity 1 — e < 10"'^ in one of the binaries. This is about 
4 to 12 times more than for triples with equivalent parameters. Orbital "flips" and collisions 
without previous tidal interaction are similarly enhanced in quadruples relative to triples. We 
argue that the frequency of evolutionary paths influenced by KL cycles is comparable in the 
triple and quadruple populations even though field quadruples are a factor of ~' 5 less frequent 
than triples. Additionally, quadruples might be a non-negligible source of triples and provide 
fundamentally new evolutionary outcomes involving close binaries, mergers, collisions, and 
associated transients, which occur without any fine tuning of parameters. Finally, we study the 
perturbations to a planetary orbit due to a distant binary and we find that the fraction of orbital 
flips is a factor of 3 to 4 higher than for the corresponding triple system given our fiducial 
parameters with impUcations for hot Jupiters and star-planet collisions. 

Key words: Binaries: close — binaries: general — planets and satellites: dynamical evolution 
and stability — stars: kinematics and dynamics 



1 INTRODUCTION 

The long-term orbital evolution of hierarchical triple syst ems of 
stars and planets is subject to Kozai-Lidov (KL) cycles (Kozail 
1 1 9621 : iLidovl 1 1 962h that transfer the angular momentum between 
the inner and outer orbits. As a result, high eccentricity and close 
pericenter passages of the inner binary can occur and physical ef- 
fects that would not be possible otherwise can significantly change 
the binary dynamics. For example, tidal friction removes orbital en- 
ergy and decreases the semi-major axis a of the inner binary, poten- 
tially explaining the exis t ence of short-period ste l lar binaries (e.g . 
Mazeh & Shaham' '1979'; 'Tokovinin et all l2006l : IWu et al.l |2007| : 
_Fabryckv & Tremaine 2007), hot Jupiters, and irregular satellites 
of planets JNesvomv et al.ll2003l) . The high eccentricities induced 
by KL cycles can also re duce the gra vitational wave merger 
timescale of compact objects JBlaes et al.l2 002 ; Miller & Hamilton 
l2002l : lThomDsoii201ll ; lAntonini & Peretsl2012l : INaoz et al.l2012ij : 
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IShappee & Thompsonll2013l : iHamers et al.ll2013l : Antognini et al., 
in preparation). 

The phenomenology of KL cycles is usually investigated at 
several levels of approximation ranging from analytic equations for 
the secular evolution of the orbital elements to fully numerical stud- 
ies. In the simplest case of the secular evolution of a test particle 
subject to a three-body Hamiltonian expanded to quadrupole order 
(hereafter TPQ limit), the KL cycles occur only if the initial mu- 
tual inclination i between the inner and outer orbits is smaller than 
cos^ iK = 3/5 (e.g. Kozai 1962). The test particle achieves a max- 
imum eccentricity of emax = \/l — (5/3) cos-^ i. In more realistic 
settings, the dynamics become more chaotic and the restrictions of 
the simplest case, especially the maximum eccentricity, do not ap- 
ply. For example, higher-order terms of the Hamiltonia n give rise to 
the eccentric Kozai mechanism iFord et al. 2000 ; Naoz et alj20Ill 
l20I3l ; lLithwick & Naozl201ll ; lKatz et alj20Ilh . which operates for 
unequal-mass inner binaries with eccentric outer orbits, giving rise 
to arbitrarily high eccen tricities and flips of the orientation of the 
inner orbit. In addition, iKatz & Dona ( 1201 2) showed with direct 
integration that for a small fraction of triples the angular momen- 
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turn of the inner binary can go from a finite value to essentially 
zero in one orbit, which can produce stellar collisions with no prior 
tidal or gravitational wave interaction. This mechanism can perhaps 
produce supemovae la by colliding two white dwarfs. Although the 
full parameter space exploration of KL cycles in triple systems with 
direct three-body integration remains to be explored, the basic prin- 
ciple that systems with lower initial | cos i \ reach higher eccentric- 
ities is still valid. 

One of the possible stable hierarchies of quadruple systems are 
two binaries on a mutual orbit. In these systems, each binary acts 
as a distant perturber inducing KL cycles on the other binary. Since 
the perturber is not a point mass, the evolution of such quadruple 
systems can in principle differ from a combination of two uncou- 
pled three-body KL processes. KL cycles in quadruples have not 
been studied previously. Recent discoveries of double close eclips- 
ing b inaries with periods very close to a 3:2 ratio (Cagas & Peicha 
uOlalKolaczkowski et alJlOlST) indicate that quadruples might ex- 
hibit new and unexplained dynamics that might be related to KL 
cycles. Moreover, a quadruple with mutual KL cycles has also been 
proposed as the origin of Tycho B as a surviving companion to the 
Tycho SN progenitor JThompson & Gould 2012; Kerzendorf et al. 
l2012h . 

In this paper, we show that for a significant part of the param- 
eter space the evolution of two binaries on a mutual orbit does not 
reduce to two independent systems consisting of a binaiy and dis- 
tant point-like perturber. As a result of more degrees of freedom 
and the mutual coupling, the two binaries evolve in concert and ex- 
perience excursions to very high eccentricities that would not be 
achieved otherwise. In Section |2] we describe our calculations. In 
Section[3] we compare the evolution of quadruples to that of triple 
stars. In Section [4] we study collisions in quadruple systems. In 
Section [5] we summarize our results and discuss the implications 
of our findings for the evolution of stars and planets. 



2 CALCULATIONS 

We modified the A'^-body code fewbody JFregeau et al.l 120040 
to simulate the evolution of quadruple stars. The integrator of 
fewb ody is based on the GNU Scientific Library dGalassi et al.l 
UOIlh . We tested Runge-Kutta Prince-Dormand (rkSpd) and the 
Bulirsch-Stoer method of Bader and Deuflhard (bsimp) and sev- 
eral conditions for adaptive timestep algorithm, all of which give 
qualitatively similar results for a few test cases. We do all of the cal- 
culations shown here with the default adaptive timestep algorithm 
with rk 8pd, which is much faster and conserves energy better than 
bsimp. We also perform the calculations using Kustaanheimo- 
Stiefel regularization JHeggid 1 1974 iMikkolal Il985h . which im- 
proves the treatment of close encounters and provides well-behaved 
long-term energy conservation. 

The dynamics of the four-body problem is inherently chaotic 
and eventually the numerical orbits deviate from the true ones with 
the same initial conditions. For calculations with absolute and rel- 
ative accuracies between 10^^^ and 10~^* we find that the system 
trajectories typically start to deviate from each other after a few 
times 10* orbits of either of the inner binaries. We continue the 
integration for longer times since we are interested in the statisti- 
cal behavior of quadruple systems. We perform all our calculations 
with absolute and relative integration accuracy of 10"^^ and we 
find that this was likely an unnecessarily high accuracy since the 
occurrence of 1 — e < 10^'' (see Section[3]l is higher only by 1 to 
3% for accuracy of 10^^^ while it is not converged for accuracy of 



10" 



. We also verified our calculations by performing simulations 



of triple stars, where the properties of KL cycles are known and we 
successfully reproduce them. Our calculations conserve total angu- 
lar momentum to about 10~^ and total energy to about 10 ~* of the 
initial value for our typical relative separations (aAB/o-A ~ 10). 
However, for large relative separations (aAB/ciA ^ 50) or very dif- 
ferent masses of the components, energy conservation is worse. We 
monitor the energy conservation of the calculations and remove all 
instances when the final relative energy conservation is worse than 
lO""". However, for our fiducial choice of accuracy, the constraint 
on energy conservation is violated by only several runs in many 
thousands and thus the statistics of our results are not affected. 

Since the parameter space of quadruple systems composed of 
binaries A and _B on a mutual orbit AB is vast, we set up the 
calculations by specifying the masses of all four bodies (tha ~ 
mi + 7712, rnB ~ raz + m^), and semi-major axes (oa, as, has) 
and eccentricities (eA, es, eAs) of the three orbits. We system- 
atically vary the angles ja and is between the angular momenta 
of orbits A and B with respect to their mutual orbit AB. The re- 
maining parameters (arguments of pericenter, longitudes of ascend- 
ing node, and mean anomalies) are chosen randomly. Typically, 
we perform 100 different random initializations for each pair of 
(cos lA, cos is). 

Due to computat ional limitations, we run the calculations for 
300 Kozai times tx l lHolman et al] 1 19971 : Urmanen et alJI 19971) of 
system A defined as 
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Higher-order effects in KL cycles become important only on 
longer timescales than ^k- For example, the secular evolution of 
the eccentric Kozai mechanism occurs on a timescale of ta/toc t 
JLithwick&NaozlbOllI : iNaoz et al.ll20Ill 1201 3l : iKatz et alj|201 ih . 
where toct is the relative strength of the octupole to quadrupole 
term in three-body Hamiltonian expansion: 



Eoct — 
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eAB 
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mi +1712 I aAB 1 — ej 



(2) 



The eccentric Kozai mechanism disappears for equal-mass binaries 
or for circular outer orbits in the secular approximation. Although 
Equations (TJ and ^ are informative for triple systems, it is unclear 
how they generalize to the quadruple case. Nonetheless, we use Ik 
to determine the length of our simulation runs and eoct serves as 
a characteristic of the triple dynamics that we use as a basis for 
understanding some quadruple phenomenology we observe. To get 
our results in finite time with the computational resources available 
to us, we limit the duration of any single calculation to 1500 s. This 
limitation becomes important only for aAB /a a ^ 20. The quadru- 
ple dynamics is different from two uncoupled three-body KL cycles 
on timescales shorter than the limits we impose on the duration of 
the calculation. Thus, our results should be very close to the "true" 
results obtained with much longer calculation times, although tech- 
nically we obtain lower bounds on the true result. We run the refer- 
ence triple calculations for the same number of f k as the quadruple 
calculations. In a number of calculations presented in this paper, we 
stop the integration when a certain value of eccentricity is reached 
or when the orbital orientation flips. 
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Figure 1. Illustration of two possible behaviors of eccentricities e and 
binary orientations cos i (measured with respect to the mutual orbit) in 
quadruple systems. The two upper panels show the two binaries evolving in 
concert and the two lower panels depict evolution where both binaries flip 
orbital orientation and go through a high-eccentricity phase. The behavior 
of both systems differs qualitatively from the expectation of two uncoupled 
triple systems. The blue lines are for system A with mi = m2 = 1 Mq, 
aj\ Ri 0.57 AU, and e^ = 0.3, while the red lines are for system B with 
ma = 1714 = lAf0, as ~ 0.35 AU, and cb = 0.1. The initial in- 
clinations are given in the plots. The mutual orbit has a/^g = 5 AU and 
EAB = 0.3. 



3 COMPARISON WITH TRIPLES 

In this Section, we investigate the properties of the quadruple dy- 
namics and we show that it is not simply a superposition of two 
independent KL cycles. Our fiducial calculation is for equal-mass 
binaries (mi — m2 = ma = 7714 = IMq) to ehminate the in- 
fluence of eccentric Kozai mechanim for the equivalent triple sys- 
tems (eoct = 0), though we consider unequal-mass binaries in Sec- 
tions |3]2] and [3]3] Summary of our quadruple calculations is given 
in Tabled 

Figure [T] illustrates two examples of behavior we observe in 
our calculations |j The two upper panels show an initial phase of 
chaotic and apparently independent evolution of the two binaries 
followed by their long-term synchronized evolution with eccentric- 
ity modulations occurring in opposite phases. Although the syn- 
chronization eventually ceases, we speculate that such mutually co- 
herent evolution is a necessary element of any the ory attempting to 
explain the near-integer period ratios in CzeV343 (ICagas & Peichal 

^ We note that the evolution of the three orbits is very close to secular in the 
sense that the standard deviations of the three semi-major axes, and angular 
momentum and eccentricity of the outer orbit around the initial values are 
< 0.5%. 



I2OI2I) and BI 108 JKolaczkowski et aljlioili . However, the cur- 
rent parameters of these observed systems dictate that such a the- 
ory must include tidal forces; we plan to investigate the interplay 
between KL cycles and tidal friction in a future work. 

The two lower panels of Figure [T] show a qualitatively differ- 
ent behavior: binary A initially shows only small eccentricity varia- 
tions compatible with its relatively large initial inclination cos ia ~ 
0.3, while binary B reaches eccentricities 1 — e_B ~ 10^^ to lO^"', 
which is expected from its initial inclination cos j_b ~ —0.05. At 
t ~ 148iK, binary B briefly flips the direction of its orbit and 
reaches 1 — es « lO^**. The evolution continues in a similar fash- 
ion until t ~ 200iK when both binaries flip their orbit orientation. 
As a result, binary A goes through many epochs of high eccentric- 
ity (1 — CA ~ 10^*) and the eccentricity oscillations in binary B 
become much shallower, reaching only 1 — e^ ~ O.L We empha- 
size that in this calculation system A had a moderate initial value 
of inclination, cosza ~ 0.3, and that a triple system analogous 
to system A with binary B replaced by a point mass would not 
achieve such high eccentricities. In the TPQ limit (mi = m2 and 
hence eoct ~ 0), the maximum eccentricity of binary A would be 
eA.max ~ \/l - (5/3 )0.3^ ~ 0.92 (Kozai 1962; I nnanen et alj 
ll997l ; lNaozetai]|2013l) . 

To characterize the quadruple properties more completely, we 
performed calculations spanning the full range of (cos ia, cos is). 
For each calculation, we recorded the highest achieved eccentric- 
ity of each binary and we show the median of these values as a 
function of (cos 1,4, cos is) using contours in Figure |2] The re- 
gion of high eccentricities, 1 — e < 10""^, extends well beyond 
cos lA ~ and cos is fa and thus a significant fraction of both 
binaries in the quadruple experience high eccentricities. In the TPQ 
limit, we would expect ^(3/5)(l -0.9992) ^ 3% ^f gj] ^^j^^_ 
ries to reach such an eccentricity. By direct integration of triple 
orbits we obtain about 7% for a.As/tSA ~ 10. This is still much 
smaller than for quadruples, where we found that 36 and 53% of 
quadruples reach high eccentricity in binaries A and B, respec- 
tively (see Section lBTl l. The region of high eccentricities is approx- 
imately bounded by the Kozai angle cosJk = ±\/3/5 of each 
binary. The asymmetry of the contours suggests that the strongest 
effect is achieved when the binaries were initially on mutually ret- 
rograde orbits. The fraction of high-eccentricity systems is higher 
for binary B than for binary A. The reason is that the perturba- 
tion to binary B from binary A and vice versa is proportional to 
a power of aA/o-AB and as/aAB, respectively, and thus with the 
same masses and uab-, the high-order perturbations induced by A 
on B are bigger than those of B on A. 

Figures [T]and|2]indicate that the parameter space of high ec- 
centricities is potentially much larger in quadruples than in triples 
due to the mutual coupling of the two binaries. Now we numeri- 
cally explore the dependence of the coupling on semi-major axes 
ratios, mass ratios, and eccentricities, and compare it to triples. 

3.1 Dependence on semi-major axes ratio 

We show in Figure[3]the fraction of quadruples with 1 — e < 10"'^ 
as a function of the ratio of semimajor axes aABJaA and aABJa-B- 
The orbital evolution time of our simulations is 300iK ~ 14 Myr 
for a A ~ 20 AU and a perturbing binary on aAB/cLA ~ 15 orbit. 
We assume that the inclinations of both binaries are distributed uni- 
formly in cos i. To directly compare the quadruples with the triples, 
we performed a series of calculations where we replaced either of 
the binaries in the quadruple with a point mass with the same total 
mass. These triple have an equal-mass inner binary and thus the ec- 
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Figure 2. KL cycles in quadraple systems composed of two binaries. At each pair of initial inclinations (cos jyi,cos is), we performed 100 randomly 
initialized calculations and recorded the highest eccentricity and the occuiTence of clean collisions. The solid lines show contours of median value of the 
highest achieved eccentricity and are labelled by logj^g(l — e). The color corresponds to the fraction of calculations that experienced a "clean" collision 
(pericenter distance at least four times smaller than any of the previous pericenters) as set by the colorbar. The blue horizontal and vertical lines mark the Kozai 
angle cos ik = ± \/3/5 for each binary. The system parameters are given in Table [T]as F6. 

Table 1. Summary of quadruple system calculations. 
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A'^ is the number of points in the (cos Ia , cos jg ) plane, which might not be distributed uniformly. Each point in (cos ja i cos is ) 
plane had 100 calculations with randomly initiaUzed arguments of pericenter, longitudes of ascending nodes, and mean anomalies. 
/a (/s) is the fraction of quadruples with binary A (binary B) reaching 1 — e < 10" 3. Symbols "— " indicate simulations where 
this quantity was not tracked. 



centric Kozai meclianism does not operate. We see from Figure [3] 
tliat for aAs/aA ~ 10 (aAB /as ~ 17) 36% (53%) of quadru- 
ples reach 1 — e < 10"'^ for binary A (B). This is by a factor 
of 4.8 (10.0) higher fraction than for otherwise equivalent triples. 
For wider separation of the binaries, the enhancement decreases: 
for aAB/dA ~ 21 (ttAB/ciB ~ 34) we find high-eccentricity 
quadruple fractions higher by a factor of 2.8 (4.9) than for equiv- 
alent triples. Wider relative separations are difficult to investigate 
numerically, but it is likely that the fraction of quadruples experi- 
encing high eccentricity asymptotically approaches the fraction of 
triples as the ratio of semimajor axes increases. Based on Figure[3] 
we estimate the high-eccentricity fraction of quadruples is signifi- 
cantly enhanced with respect to triples for aAB/dA ^ 50. 



3.2 Dependence on mass ratios 

So far we have discussed only quadruples composed of equal- 
mass stars. The strength of the interaction depends also on qa = 
m2/mi, qB ~ mi/m-j,, and qab ~ mA/rriB- We expect that the 
mutual interaction vanishes for two stars orbited by a test particle 



(qA — >■ 0, gs — >■ 0), when the dynamics reduces to two indepen- 
dent KL cycles. To see how non-equal mass quadruples behave we 
ran a calculation with qA = 3/4, gs = 1/2, qAB = 7/3 (entry 
B6 in Table[TJ. As expected, the fraction of calculations where sys- 
tem A reached high eccentricity is lower (20%) than in the fiducial 
case because of the lower mass of the perturbing binary niB ■ Inter- 
estingly, the fraction increases for binary B (77%) and essentially 
all systems with cos^ is 5C cos^ zk = 3/5 reach very high ec- 
centricity. The reason is again that aA/o-AB is significantly larger 
than aB/aAB, and also tua > m.B. In triple systems with different 
inner binary masses, eccentric Kozai mechanism becomes impor- 
tant (foct > 0) and induces high ecc entricities over a wider range 
of cos i than in e qual-mass systems (Fordetal. 2000!: iNaoz et al.l 
|201 lll2013l : lLithwick & Naoz 201 1 : Katz et al. 201 1 ). We verified 
by direct integration of triple systems that for our unequal -mass 
system, the fraction of quadruples reaching 1 — e < 10~^ is a fac- 
tor of about 3.0 (binary A) to 8.2 (binary B) higher than for corre- 
sponding triple systems defined by replacing one of the binaries by 
a point mass. 
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Figure 3. Fraction of systems with 1 — e < 10""^ as a function of the 
semimajor axis ratio showing that more quadruples than triples reach high 
eccentricity. The circles connected by a line indicate results for the two bi- 
naries A (solid blue) and B (open red) of quadruple systems. The tiiangles 
and diamonds connected by a line indicate results from systems where ei- 
ther binary Aor B was replaced by a point mass with the same total mass. 
The fiducial system parameters are given in Table [T]in lines F6 and F12. 
Three overlapping circles labelled with grey arrows denote binaries with 
a A = a-B (entries A6eO, A12e0, and A24e0 in Table[T) The grey horizontal 
dashed line indicates fraction of binaries with cos^ i ^ cos^ jk assuming 
uniform distribution of cos i. 



3.3 Planets and orbital flips 

Systems of a star and planet on a mutual orbit around a distant bi- 
nary are an important example of quadruples with an extreme mass 
ratio. The planetary orbit can experience the eccentric Kozai mech- 
anism leading to high eccentricity and orbital flips, which can po- 
tentially explain the abundance of observed retrograde hot Jupiter s 



JNaozetaLkOUl , |2013| ; iLithwickfe Nao3l201ll : Ixatz et alj|201lh . 
In Figure|4] we show the flip fraction of planet orbits evaluated 
as a change of sign of the projection of the star plus planet angular 
momentum to the total angular momentum of the quadruple. For 
the parameters of this particular system and assuming a uniform 
distribution for cos i of both binaries, we find that about 66% of 
all planetary orbits flip their orientation. How does this compare 
to triple systems where the distant binary is replaced by a single 
star? The flip fraction of triples presented in the upper panel of Fig- 
ure |4] shows clear signs of the eccentric Kozai mechanism as the 
broad component of the peak at —0.5 < cosia ^ 0.3. Nonethe- 
less, the total flip fraction of ~ 21% is still ~ 3.2 times lower 
than for the quadruple system. The eccentric Kozai mechanism op- 
erates on a longer timescale, iK/toct ~ 40tK, for our choice of 
dA/dAB ~ 0.072 and cab = 0.3 and we thus ran our triple calcu- 
lation for 1500iK ~ 38fK/eoct. For our fiducial calculation length 
of 300tK ~ StK/^oct, the flip fraction was about 2% lower, sug- 
gesting that our comparison of quadruples and triples is robust[j 

^ The triple flip fraction in Figure |4] does not look completely converged 
despite 100 random initializations at each cos i^. It is possible that a sig- 
nificantly higher fraction of triples with — 0.5 < cos iA 5: 0-3 would even- 
tually flip their orbits given integration times that are orders of magnitude 
longer than what we are able to achieve. The flip fraction would be ^ 40% 
if all such orbits flipped, but this is still 1.7 times less than what we obtain 
for quadruples. It is also quite likely that the quadruple flip fraction would 
increase as well for very long integration times. 
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Figure 4. Orbital flip fractions for star plus planet systems orbiting another 
binary. The lower panel shows the flip fraction as a function of the incli- 
nation of the star plus planet (cos ia) and star plus star (cos is) binaries. 
The upper panel shows the flip fraction integrated over all orientations of 
the star+star binary (black line) in comparison to a tiiple system, which has 
the outer binary replaced by a point mass with the same total mass (green 
line). The system parameters are given in Table[T]as sjss. 



Although the system in Figure|4]is not representative of typical hot 
Jupiter systems (e.g. Wu et al. 2007; Naoz et al. 2012b) and a full 
exploration of the parameter space is beyond the scope of this pa- 
per, our results suggest that it is much easier to place a planet on a 
close and possibly retrograde orbit around a star that orbits a rela- 
tively distant stellar binary. Similarly, star-planet collisions are also 
more likely in quadruple systems than in triples. 



3.4 Dependence on eccentricities 

Finally, the coupling between the two binaries depends also on 
the remaining orbital parameters that we either specified explic- 
itly as ba, e_B, eAB, or marginalized over by randomly initializ- 
ing the calculations (arguments of pericenter, longitudes of ascend- 
ing node, mean anomalies). In order to characterize the effect of 
the initial eccentricities of the orbits, we performed several calcu- 
lations with various choices of eccentricities and we did not find 
any significant changes (for example A6 and A6e0 in Table [TJ, 
except that high initial cab gives a higher fraction of quadruples 
that become dynamically unstable. To get a cleaner view of the ef- 
fect of eccentricities and to test our code, we performed several 
calculations where all the explicitly specified parameters of the bi- 
naries A (blue solid circle) and B (red open circle) were set to 
identical values (entries A6e0, A12e0, and A24e0 in Table [TJ and 
we show the results with overlapping circles in Figure |3] Interest- 
ingly, the fraction of quadruples composed of two identical bina- 
ries that reach 1 — e < 10^^ starts decreasing considerably only 
for aAB/ciA ^ 30, which is later than for the fiducial case we 
investigated previously. Due to numerical difficulties with energy 
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conservation at large relative separations we did not investigate this 
potentially interesting issue further. 



4 COLLISIONS 

The KL cycles in triples can change the angular momentum of the 
inner binary in such a way that it goes from a finite value to es- 
sentially zero in a single orbit. As a result, stars that previously 
were too far from each other to interact through tides or gravita- 
tional wave emission will collide (Katz & Dong 2012). However, 
for equal-mass binaries this occurs only in a few percent of sys- 
tems where the perturbing body initially orbits almost perpendicu- 
larly with respect to the inner binary (cos i ^ 0). In Section[3] we 
showed that quadruples reach high eccentricity for much larger part 
of the parameter space than triples. We now investigate whether the 
same is true for collisions. 

In a subset of our calc ulations, we implemen ted a collision de- 
tection algorithm based on iKatz & Dong (2013). Since our calcu- 
lations are dimensionless, we assume that a clean collision without 
any previous tidal or gravitational wave interaction occurs when a 
pericenter distance of the binary is at least by a factor of 4 smaller 
than any of the previous ones. The eccentricity at this pericenter 
sets the radius of the stars relative to the semi-major axis. The up- 
per limit on stellar radius for such "clean" collisions to occur is set 
by the time for quadruple KL cycles to develop, which scales with 
fK, relative to the lifetime of the stars or the Hubble time. 

Color pixels in Figure |2]show the fraction of binaries that ex- 
perienced a clean collision for our fiducial calculation. Similar to 
the cases of eccentricity and orbital flips discussed in Section [3] 
the clean collision fraction is significant over a large chunk of the 
parameter space and does not require any fine tuning of the incli- 
nations. In our fiducial case displayed in Figure|2] about 19% and 
25% of binaries A and B, respectively, experience clean collisions. 
This is much higher than for corresponding triples, where we find 
4% an d 2%, quite close to the values reported by iKatz & Dong 
( 120121) for a similar ratio of semi-major axes. The distribution of 
1 — e at the moment of clean collisions starts at about 10"'' and 
peaks at about 4 x 10"'', which permits clean collisions of main- 
sequence and smaller stars in a Hubble time, as found also by 
iKatz & Don3 12012|) . 



5 DISCUSSION AND CONCLUSIONS 

Our numerical calculations of orbital evolution of quadruple sys- 
tems shows that their KL cycles are qualitatively and quantitatively 
different from triples. Although the binarity of the KL perturber 
would seem to be a higher-order effect, it significantly increases the 
fraction of close passages of stars in quadruple systems compared 
to equivalent triples. We now explore some of the implications of 
our findings. 

We showed that high eccentricities (1 — e < 10~^), orbital 
flips and stellar collisions occur with a frequency ~ 3 to 12 times 
higher in quadruples than in triples for aAs/aA ~ 10 to 20 
and we estimated that the frequency is still significantly higher at 
dAB/dA ~ 50. The strength of quadruple KL cycles in both bina- 
ries is maximized for equal masses. However, quadruples are also 
intrinsically less frequent in the field than triples. Thus, to assess 
the absolute occurrence of high eccentricities in the field population 
of quadruples and to compare it to triples, we need to know their 
relative frequencies, and the distributions of semi-major axes and 



mass ratios. The best available source is the catalog of lTokovininI 
1 1 19971. 2008). By analyzing the distance distribution of triples and 
quadruples composed of two binaries, Tokovinin ( 2008) found that 
such quadruples are by a factor of ~ 5 less frequent t han triples 
in the field. The median of aAs/aA (ciAB/aB) in the iTokovininI 
(2008) sample is 100 (150) and ~ 40% (- 30%) of quadruples 
have aAs/fflA < 50 (aAs/as < 50). For triples, the median semi- 
major axis ratio is --^ 80 and '^ 45% of systems have a ratio smaller 
than 50. The medians o f qa, <1b and qAB of 81 quadruple systems 
in the lTokovininI ( 1200 81 ) sample with mass estimates are 0.64, 0.70, 
and 0.70, respectively. To summarize, the field sample of quadm- 
ples appears to be composed of stars with similar masses and the 
relative distribution of semimajor axis ratios in quadruples is no t 
significantly different from triples. Although the lTokovininI ( 120080 
catalog is subject to many selection effects, the lower intrinsic fre- 
quency of quadruples in the field is approximately compensated for 
by the higher efficacy of their KL cycles. Quadruples are thus an 
order unity correctio n to any population synthesis result involving 
KL cycles in triples JHamers et alj2013h . 

In principle, KL cycles in quadruples will act on both binaries 
while in a triple only the inner binary will be affected; a quadruple 
thus potentially produces twice as many "interesting" outcomes. In 
reality, after the first clean collision, merger, or close encounter in 
one of the binaries of the quadruple system, this binary will become 
a single star or its semimajor axis will significantly shrink due to 
tidal dissipation and the subsequent evolution will then proceed ef- 
fectively as a triple star, with the usual properties of KL cycles in 
triples. Assuming that the initial distribution of the ratio of semi- 
major axis of the inner and mutual orbits was the same in quadru- 
ples and triples, the higher efficiency of KL cycles in quadruples 
evolution will lead to a relative lack of quadruples with small semi- 
major axes ratio. T his is supported b y higher median of aABJaA 
and aAB /iB in the lTokovininI ( 120081) sample, however a two-way 
Kolmogorov-Smimov test on the quadruple and triple semimajor 
axis ratio distributions does not reveal any statistically significant 
differences. 

The lower frequency of quadruples relative to triples in the 
field might also be an effect of highly efficient quadruple KL 
cycles: the fraction of mul tiple stars in star-forming regions is 
higher than in the field ( e.g.lLeinert et al. 1 19931; lOhezetal .ill 9971; 
iKohler & LeinertI 1 19981 ; IChen et alj 120131) and the frequency of 
triple and quadruple stars might be comparable (e.g. ICorreia et al.l 
I2OO6'). Thus, quadruples may be a significant source of triples. 
Furthermore, stellar evolution and associated mass loss will lead 
to expansion of orbits that will be more pronounced for the in- 
ner binaries than for the mutual orbit and thus yield smaller 
dAB/oLA and aAB/cLB- On the other hand, when the mass ra- 
tio of the binaries evolves farther from unity the quadruple KL 
cycles become weaker (Section |3. 2) . Depending on how exactly 
the mass loss occurs, quadruple systems that did not experi- 
ence KL cycles due to large separation of the binaries and un- 
favorable orientation of the orbit for the usual triple KL mech- 
anism might start experienci ng strong KL cycles as they evolve 
dShappee & ThompsonluOlSh . Additionally, a supernova or a neu- 
tron star kick in one member of the quadruple may also lead to in- 
teresting results dShappee & Thompson l2013l) . Finally, interesting 
evolution in triple systems of stars and planets occurs also with K L 
cycles not present dPerets & Kratteil2012tlKratter & Peretsll2012h . 

The evolution of quadruples can produce outcomes that would 
not be achieved by triples and can involve main sequence stars, 
white dwarfs, neutron stars, and their combinations. Some out- 
comes include pairs of close binaries or blue stragglers on a wide 
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orbit dPerets & Fabrvckvl l2009f) , blue stragglers orbiting a close 
binary, massive white dwarfs orbiting a (close) binary or a blue 
straggler, etc. Collisions and mergers of stars produce transients 
that can be observed over great distances . Potential examples in- 
clude V838 Mon and V1309 Sco (e.g. 'Soker & Tylenda' '2003!; 
Tvlenda & So ke] [20061: b lebbeek et al. 2008; Glebbeek & Pols 
20081 ; ISmitHl201 ll ; iTvlenda et al.ll201 Ih . Additionally, supemovae 



la c an occur as a re s ult of a collision o f two white dwarfs 
(e.g. iBenz et"d] Il989l; iRosswog et alj |2009|; JRaskin et alj l2009l 



2010l;lLoren-Aguilar et al.l2010l;lHawlev et al.l2012l;lKushnir et al.l 



20131) . The KL mechanism allows such collisions to happen in the 
field, not only in dense stellar environments, but only a few per- 
cent of triples can collide due to severe restriction s on the incli- 
nation of the perturbing body (IKatz & Dongll2012r) . However, we 
find that quadruples produce collisions generically for a large frac- 
tion of the parameter space and thus do not require fine tuning 
in the orientation of the orbits. Furthermore, the quadruple evo- 
lution allows for peculiar stars such as blue stragglers, close bi- 
naries, or massive white dwarfs to remain at the explosion site. 
Specifically, our results on KL cycles in quadruples increase by 
a factor of few the probability of having a metal-poor A star (blue 
straggler) in the Tycho sup ernova remnant (IKerzendorf et al .120 121 ; 
[Thompson & Goulq|2012i) . The higher efficiency of quadruple KL 
cycles might give important constraints on the future evolution 
of the kn own quadruple s ystems such as t hose i nvestigated by 
iTokovinin et al. ( 2003) and iHarmanec et alj ( 120071) . The new dy- 
namics of quadruple systems together with tidal dissipation might 
explain the 3:2 period ratios of some double e clipsing binaries 
dCagas & Peichall2012l ; lKolaczkowski et alj|2013l ), but more work 
is required to understand these systems. 

Finally, our results imply that stars orbiting a distant stel- 
lar binary have significantly higher probability of secular changes 
to their planetary orbits (Fis^|4) than in systems where the dis- 
tant body is a point mass (^ u & Murray 2003; Wu et al. 2 0071: 
Fabryckv & Tremain e 2007; Naoz et al.l 1201 iL 120131; iKatz etaU 



pry' 

m 



201 IclLithwick & Naoz..201 1.) . Without a more detailed investiga- 



tion of the parameter space, it is not clear whether the enhanced 
KL cycles in such systems would lead to a higher fraction of hot 
Jupit ers (possibly on re trograde orbits) or more star-planet colli- 
sions. [RoeUetalJ J2012i) list nine exoplanet host stars that are mem- 
bers of triples that are in the hierarchy considered here. For most 
of these systems, the separation of the host star from the distant 
binary relative to its separation is within the range considered in 
this paper, aAB/aB ^ 50. Planets in these systems are not hot 
Jupiters, but the secular perturbation timescale would be large due 
to the size of the stellar orbits, which are often visually resolved. A 
planet subject to quadruple KL cycles on a short timescale would 
be the progenitor of the hot Jupiter orbiting HD 188753 A, member 
of a triple stellar system (Konacki 2005). However, thi s planet was 
not independently confirmed dEggenberger et al.ll2007h . 
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